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ON INVARIANT MEASURES OF STOCHASTIC 
RECURSIONS IN A CRITICAL CASE 

By Dariusz Buraczewski^ 
University of Wroclaw 

We consider an autoregressive model on R defined by the recur- 
rence equation X„ — AnX„-i + B„, where {{Bn, A„)} are i.i.d. ran- 
dom variables valued in R x R"'" and E[logyli] — (critical case). It 
was proved by Babillot, Bougerol and Elie that there exists a unique 
invariant Radon measure of the process {Xn}. The aim of the paper 
is to investigate its behavior at infinity. We describe also stationary 
measures of two other stochastic recursions, including one arising in 
queuing theory. 

1. Introduction. We consider the following random process on M: 

where the random pairs {Bn,An) in M x M'^ are independent, identically 
distributed (i.i.d.) according to a given probability measure /x. This process 
is called sometimes a first order random coefficients autoregressive model. 
It appears in various applications, especially in economy and biology; see, 
for instance [1, 18, 25] and the comprehensive bibliography there. 

Many properties of the random process {-'^n} have been studied under 
the assumption 

(1.1) E[log^i]<0. 

Then, there exists a unique invariant probability measure v of the process 
{Xn} ([17]; see also [12]). This is a measure such that 

^ *!/(/) = !/(/), 
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for any positive measurable function / on M, where 

/^*^(/)= / / f{ax + b)diy{x)dfi{a,b). 
JMxR+ Jm. 

Under some further assumptions, the tail behavior of v has been described 
by Kesten [17]. He has proved that 

(1.2) >x})~Cx"" asx^+oo 

for some positive constants a and C. Kesten's proof was later essentially 
simplified by Grincevicius [15] and Goldie [11]; see also [5, 8, 13, 16, 18, 19] 
for related results. 

In this paper we are going to study the so called "critical case," that is, 
the case when 

(1.3) E[logAi]=0. 

Then there is no finite invariant measure. However, it has been proved by 
Babillot, Bougerol and Elie [1] (see also [3, 4]) that there exists a unique 
(up to a constant factor) stationary Radon measure of the process {Xn}- 
Moreover, they have described the behavior of v at infinity, proving that, 
for any given positive a and f3, 

^{{ax, fix]) log{P/a) ■ L^{\x\) as x ^ ±oo, 

where are slowly varying functions. 

The aim of this paper is to show that under the hypothesis that the 
measure fi is spread out and has some moments, slowly varying functions 
are constant. Indeed, we show 

(1.4) iy{{ax,(3x]) ^log{P/a) ■ C± as x — > ±oo 
(Theorem 2.2). 

If the measure is of a very specific form, that is, fi = jit, where {nt}t>o 
is a one-parameter semi-group of probability measures whose infinitesimal 
generator is a Hormander type differential operator on R"' X M"*" (or, more 
generally, on solvable groups of NA type), then the invariant measure u has 
smooth density m and in the situation corresponding to the critical case, 

(1.5) x'^m{x ■ t) c{t) as X — > +00, 

for any t ^W^\ {0}, where x • t is an appropriate dilation, [6, 7]. Of course, 
for these particular measures /x, the present result (1.4) follows from (1.5). 

In the contractive case (1.1) it was observed by Goldie [11] (see also 
Grey [13]) that the same problems can be investigated, if the random linear 
transformations t^At + B are replaced by a general family of transforma- 
tions {^(t)}j>o, provided that, for large values t, ^(t) is comparable with 
At. Goldie studied several random processes, proved existence of stationary 



STOCHASTIC RECURSIONS IN CRITICAL CASE 



3 



probability measures and described their behavior at infinity that coincided 
with (1.2). In the critical case we examine a model due to Letac [21]: 

X'^ = Bn + An max{X^„i , C„,}, 

where (An, Bn,Cn) € IR"*" x M x M'^ are i.i.d. Then, existence and unique- 
ness of a stationary measure follows from the theory of locally contractive 
stochastic dynamical systems due to Benda [3]. We shall prove that also in 
this case the tail behaves regularly and satisfies (1.4) (Theorem 5.1). In par- 
ticular, the result holds when = and we consider the following random 
process: 

<' = max{^„<'_i,A.}, 

where {An,Dn) € M"*" x M"*" are i.i.d. The process {^^} is called the extremal 
random process and it plays an important role in modeling of the waiting 
time for a single server queue [10]. 

The structure of the paper is as follows. In Section 2 we describe the 
autoregressive process in the critical case and state our main result. Theorem 
2.2. Next, in Section 3 we describe results concerning solutions of the Poisson 
equation and their asymptotic behavior. In Section 4 we conclude the proof 
of Theorem 2.2. Finally, in Section 5 we investigate the Letac's model. 

2. Random difference equation Xn = AnXn—i -\- Bn. 

2.1. Main theorem. Given a probability measure ^ on M x M"*", consider 
the following Markov chain on R: 

Xo = 0, 

(2.1) 

Xn — AnXn~\ -\- Bn, 

where {Bn, An) are i.i.d. random variables with values in M x R+, distributed 
according to /i. 

The process {Xn} is best defined in the group language. Let G be the 
"ax -|- 6" group, that is, G = M x M+, multiplication being defined by 

{b,a) ■ {b' ,a') = {h + ah' ,aa'). 

G acts on M by 

{b,a) o X = ax + h, {b,a)^G, x G M. 

We sample {Bn,An) G G independently according to a measure /i and we 
write 



Xn = {Bn,An) o = (B„, A„) • • • {Bi,Ai) o 0. 
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This paper is about the critical case, that is, our main assumption is 

Elogyl = 0. 

We are interested in the asymptotic behavior of the (unique) invariant Radon 
measure of the process that is, the measure on R satisfying 

(2.2) ;,*gK/) = K/), 
for any positive measurable function /. Here 

/^*Gi^(/)=/ / f{ax + b)du{x)dfi{a,b). 

JrxR+ Jr 

(We write *g for the convolution induced by the action of G on M defined 
as above, in order to distinguish it with the convolution on M, i.e., denoted 
by *R.) 

Existence and uniqueness of such a measure i' is due to Babillot, Bougerol 
and Elie [1] (see also [3, 4] for some comments), who proved the following 
result: 

Theorem 2.1 ([1]). Assume 

(2.3) ElogA = 0; 

(2.4) A^l; 

(2.5) F[Ax + B = x]<l for all x£ M; 

(2.6) E[(|logA| +log+ |S|)2+'5] < oo for some 6 >0. 

Then there exists a unique (up to a constant factor) invariant Radon measure 
u onM of the process {Xn}- Moreover, if 

the closed semigroup generated by 
the support of pL is the whole group G, 

then there exist two slowly varying functions and on M"*" such that, 
for any a, /3 > 0, 

(2.8) v{{ax,l3x])r^\og{l3/a)-L^{\x\) as x^±oo. 
In particular, 

(2.9) / {l + \x\)"^ dv{x) <oo, 

Jr 

for any 7 > 0. 
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Our aim is to study more precisely the behavior of v at infinity. We shall 
prove, under additional assumptions, that the functions and L" are just 
constants. 

Define a probability measure [xa on M"^ being the projection of onto the 
second coordinate [ia = '^2{fJ'), that is, for any Borel set U contained in M^, 
we put ha{U) = fi{M X U). Recall that a measure on M"*" is called spread- 
out if for some n, its nth convolution power has a nonsingular component 
relative to the Haar measure on M+. 

Our main result is the following: 

Theorem 2.2. Suppose that assumptions (2.3)-(2.7) are satisfied and, 
moreover, there exists a positive constant 6 such that 

(2.10) EA^<oo and EA'^ < oo; 

(2.11) E\B\^<oo; 

(2.12) fiA is spread-out. 

Then, for any positive numbers a, (3, such that a < (3, we have 
lim ^{{ax, (3x]) =log{f3/a) ■ C-^-, 

X— >+oo 

(2.13) 

lim i^((/3x, ax]) = log(/3/a) • C_. 

X— » — oo 

2.2. Sketch of the proof. Fix two positive numbers a < f3, and define a 
function / on M: 

(2.14) /(x) = z.(ae^/3e^. 

Define a measure /2 on M: 

fl{U) = fiA{{x : - logx G [/}), 

for any Borel set U. By (2.3), the mean of p, is equal to 0. 
Define convolution of a function g with a measure on M: 



r]*Rgix)= / g{x + y)r]{dy). 
Jr 

The key observation is the following lemma. 

Lemma 2.3. The function f satisfies the Poisson equation 
(2.15) i2*^f{x) = fix) + i;{x), 

where 



(2.16) ip{x) 



ae^ I3e^ 
a ' a 



( ae^ — b Pe^ — b 
\ a a 



dfj,{b, a). 
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Proof. We have 




= / ; dfi{b,a) 



and, hence, 




ae^ — b (5e^ — b' 



dfi{b, a) 



JRxM+ 



a 



a 



r r fae- f3en 

= f{x)+ / z/ ; 

JRxR+ L \ a a _ 



— V 



( 



ae^ — b Pe^ — 61" 



a 



a 



dii{b, a) 



= f{x)+ip{x), 



which proves (2.16). □ 

The Poisson equation on M was studied in the 1960s (see next section for 
more comments). It is weh known that if the function -0 is good enough, 
there exists a formula describing solutions of (2.15). This is not our case, 
however, in the next section we prove that even if ip does not satisfy classical 
hypotheses, under some other assumptions (Theorem 3.3), one can obtain 
enough information to describe the behavior of the function / at infinity. In 
Section 4 we prove that tp possesses all the required properties and conclude 
the proof of Theorem 2.2. 

3. Poisson equation. 

3.1. Classical results. Given a probability measure on R that is cen- 
tered, that is, J^xdfi{x) = 0, and a function cj) on R, consider the following 
Poisson equation: 

(3.1) fl*u f{x) = f{x) + <j){x) dx a.s. 

Positive solutions of this equation were studied in a number of papers (see 
[2, 22, 24, 26]). For our purpose, we use the result of Port and Stone [24], 
who have considered the situation when the function (p is bounded and 
compactly supported. Then, assuming that the measure p, is spread out, 
they have found an explicit formula for solutions of (3.1) that are bounded 
from below. Weaker assumptions on fl were imposed by Ornstein [22] and 
Baldi [2]; however, they studied only positive functions (p. 
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Port and Stone [24] define a potential kernel A, which can be explicitly 
computed for (p G ^ is the class of bounded measurable, compactly sup- 
ported functions). The kernel is of the form 

(3.2) A(l){x) = a*R(l){x)- H2*R<l){x)+bJ{(p)-(f){x), (/) G $ 

([24], Theorem 7.1), where is a finite positive measure, b is an appro- 
priately chosen constant, J{(j)) = 4>{x) dx, a is a continuous function [23] 
satisfying 

(3.3) lim {a{x -y) - a{x)) = ^a~^y, 

x—>±oo 

where = J^x'^j2{dx), and moreover, the convergence is uniform w.r.t. y in 
compact sets. 

The potential A provides solutions of the Poisson equation: 

Theorem 3.1 (Port and Stone [24], Theorem 10.3). Assume that the 
probability measure p, is centered, spread-out and its second moment is finite. 
Then, if (p is an element of^, all solutions of Poisson equation (3.1), which 
are bounded from below, are of the form 

(3.4) f(x)=A(P(x) + ^^^^x + d, dx a.s., 
where d is any constant and |c| < 1. 

We observe that (3.2) and (3.3) imply that for large x and |c| < 1 the 
absolute value of the term '^'^^f^ x is dominated by a *k (j), so f{x) defined by 
the right-hand side of (3.4) is bounded from below. 

The Poisson equation can be solved for a more general class of functions 
and all solutions are given by a formula that coincides with (3.4). Indeed, 
the following result holds. 

Theorem 3.2. Assume that a function (f) is bounded, continuous and 
satisfies f^ \4){x)x \ dx < oo, then A(p is a well defined continuous function 
and all solutions of the Poisson equation (3.1) bounded from below are of 
the form (3.4). 

The proof of the foregoing theorem follows the classical path; we present 
the details of the argument for the reader's convenience in the Appendix. 

3.2. Behavior at infinity. Now, we are going to present a result that 
will be the main step in the proof of Theorem 2.2 and will also be used in 
our study of the asymptotic behavior of the invariant measure v for other 
recursions considered in the paper. 
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Let be a Radon measure on M. Assume that for every positive 7 there 
exists a constant C such that 

(3.5) i^{0,x]<C{l + x^), x>0. 

Let fa,p{x) = v{ae^ , I3e^]. Suppose fa^p satisfies the Poisson equation 

(3.6) /i*R/a,/3(a;) =/a,/3(x) +^a,/3(x), X € M, 

where the measure Jl satisfies the assumptions of Theorem 3.1 and also 

POO 

(3.7) / e^^ djl{x) < 00 



for some 7 > 0. 

Theorem 3.3. Assume 

/•CO 

(3.8) / \ipa,i3{x)^\dx < 00; 



00 

1 



(3.9) J \ipa,i3{x)\ dx < QO. 

We let 

Cl,l3 = - {x + l)V'a,/3(2;) dx. 



Suppose that 

(3.10) the functions {a, (3) 1— ^ C^^ ^ and (a,/3) i— >C^^ are continuous. 
If a function f^^/s satisfies the Poisson equation (3.6), then 



Ct2 



hm 

x—>+c 

Moreover, if p = 0, then 

2(^2 

lim fa,/3{x) = j^. 

a;— >+oo (J 

Our aim is to solve the Poisson equation (3.6). Assumptions (3.8)-(3.10) 
are not sufficient to do it directly; therefore, we are going to use a smoothing 
operator to obtain a "smoothed" version of (3.6). This will satisfy all the 
hypotheses of Theorem 3.2. We prove that the solution of the regularized 
equation behaves regularly at infinity, and finally, we obtain the description 
of the asymptotic behavior of the function fa,i3- 
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Consider the function K{t) = e^*l[o,oo)(*) define a smoothing opera- 
tor (cf. [11, 15]): 

(3.11) g{t) = K*^g{t)= f e-^'-''^g{u)du. 

Lemma 3.4. The functions ipa,f3, fa,p o,nd {fi *m. fa,i3y o'^e well-defined. 
Moreover, {p. *r fa,i3y= P- *R fa,i3- 

Proof. By (3.8), Va,/? is integrable, hence, ipa,i3 as the convolution of 
K with ipa,i3 is well defined. Next, for every t, we have 

t ft 



faAt)= e~(*-")i^(ae";/3e"]d'u<zy(0;/3e*] / e'^*"") dn < oo. 

J — OO J — OQ 

Observe, that by (3.5) and the Tonelli theorem, we obtain 
(/i *R /a,/3)t*) = / e~(*~") / f^Au + x)dJi{x)du 



oo 

t+x 
-oo 



e-^*+"-"Va,/j(^)d«d/i(^) 
< / z^(0,/3e*+^](i/i(x) <C / (l + e^(*+^))(i/i(x). 



which, in view of (3.7), is finite. Finally, the Fubini theorem implies the last 
conclusion of the lemma. □ 

Applying the smoothing operator to both sides of (3.6), we obtain a new 
Poisson equation: 

(3.12) A*M/a,/3(a^) =/a,/3(a^) +V'a,/3(2;), X G M. 

Moreover, in this case the function -00,^^ is good enough to yield a description 
of solutions of the equation. 

Lemma 3.5. The function tpa,f3 satisfies assumptions of Theorem 3.2, 
and moreover, it vanishes at infinity. 

Proof. Of course, the function is continuous. To prove bounded- 
ness of V'o,/?, it is enough to prove that its limits at +oo and — oo exist. In 
fact, as we are going to prove, the limits are equal zero. We have 

lim |V'Q:,/3(t)| < lim / |V'q,/3(u)| = 0, 

t— »— oo I— >— oo 
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oo 



lim |^c.,/3(t)| < lim / e ")l[o,oo)(i - w)|?/'a,/3(w)| o!?^ 



oo 



lim e ""hiQ ){t-u)\tpa,f3{u)\du = 0. 
The validity of the equality follows from Lebesgue's theorem. We also have 

roo rx 



-oo J — oo 
oo 



(3.13) 

rO 



oo 



\'4'a,i3{x)x\dx < I I Ixe*^* ^'^'tpa,i3{t)\dtdx 

e*|V'a,/3(*)l • {^J Ixle'^'dxj dt 
e*|^a,/3(t)|(2-(t + l)e-*)(it 

oo 
/■oo 

+ / (t + l)|Va,/3(t)|dt. 





The value above is finite, by (3.8) and (3.9). □ 
Proposition 3.6. We have 

(3.14) lim kf^ = %L_ 

X— »+oo X a 

Moreover, if p = 0, then 

(3.15) hin f\p{x) = ^^. 

Proof. The function /q,^^ is positive and satisfies the Poisson equation 
(3.12), moreover, ipa^fj satisfies the assumptions of Theorem 3.2 (Lemma 
3.5), therefore, Atp^^fs is a well-defined continuous function and 

(3.16) faAx) = A-ii^aAx) + '^'^^'^"'^'^ x + d, dx a.e. 

Notice that both sides of the foregoing equation are continuous functions, 
therefore, it is satisfied for all x € M. Next, observe that because is a Radon 
measure, for any positive e, there exists a negative M such that fa,p{u) < s 
for u < M . Therefore, for x < M, 



\faA^)\< e-(^-")/,,;3(^/)dn<e, 

which implies 

lim faij{x)=0. 
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(lix) 1 

lim = ±— ^ and \a{x — y) — a{x)\ < C{\y\ + 1) 



for some constant C and every x € M. Therefore, applying (3.16) and Lemma 
3.5, we obtain 



0= lim fa,i3{x)= lim 



■d 



lim 

X— > — oo 



+ d 



bJ{lpa,l3) + d 



+ lim 

x—*—oo 



{a{x -y)- a{x))'il^aAy) dy + ( a{x) + -^x ) J(^a,/3) 



- bJilpa,^) + d + ^ '4)a,f3{y)y dy + J{ll)a,p) ■ 1™ X 
a .IV! x~*~ca 



a(x) c 



a' 



which implies J{ipa,f3){c — 1) = and 
In the same way we compute 



lim 4^ = J(^„^^). lim 



a{x) 



2J(4,/3) 



x^+co X x^+oo\ X 

but one can easily prove J{il'a,i3) = •^(V'a,/?) = C^pi which gives (3.14). If 
<i 



Cl . = 0. then 



lim fa^fiix) = 2 

X—f+OO (J 



4^aAy)ydy, 



and observing 

/ '(PaAy)ydy= / {y + ^)i^aAy)dy, 

we obtain (3.15), which completes the proof. □ 
Proof of Theorem 3.3. We have just proved 

hm (j;logx) / fa,i3 {log t) dt = 

x—>+co Jq a 
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Fix e > and consider 1 < 5 < 1 + e. We have 

[S-l) '——^ >(2;logx / fa,/3(logt) dt ^ — ((5-1). 

Therefore, 

(3.18) ^■^■^^ /.,(i^.)Mlog^) ^gg^. 

x^+oo log X 

Analogously, taking < 5 < 1 , we prove 

(3.19) limsup:^-Ml±f)l!^<!^. 
Comparing (3.18) and (3.19), we obtain 



< lim inf — — < lim sup — < 



for every e > 0. Finally, passing in the foregoing inequality with e to zero 
and applying (3.10), we conclude 

lim t^J^^ 

rr^+oo X a 

which completes the proof. The same argument justifies the second part of 
the theorem. □ 

4. Proof of Theorem 2.2. In order to prove Theorem 2.2, it is enough 
to check that all the assumptions of Theorem 3.3 are satisfied and compute 
explicitly constants. We begin with a simple lemma. 

Lemma 4.1. For any 7 < 1, there exists xq € (0, 1) such that the measure 
I' = ^xo *R satisfies 

(4.1) [' J-di)is)<oo. 

J~i \s\' 

Proof. We have 

rl 1 /•! /•! 1 

■dv{s)dx= / / 1 —dxdv{s) 



Jo \s-x\i Jo Jo \s-xp 

1 rl 



<[ [ j^dxdu(s) = C ■i^[0,l]<oo, 
Jo J-i Ixu 



lo J-i \xp 
therefore, there exists xq G (0, 1) such that 

(4.2) L^du{s)<^, 

Jo \s-xoP 
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which imphes the result. □ 



Define the following measure on G: fl = S(^xo,i) *G *G '5(-xo,i)) that is, 
= ;^((— xo,l) • U ■ (xq,!)), for Borel sets U CG. Then fl satisfies all 
the assumptions of Theorem 2.2. Hence, there exists a Radon measure on M 
invariant under fl and one can easily check that this is exactly the measure u 
(up to a constant factor) defined in the lemma above. Of course, the behavior 
of both measures v and D at infinity is the same. Therefore, it is enough to 
prove Theorem 2.2 for i). From now we shall consider the measures p, and 
i> instead of /x and i'. However, to simplify our notation, we just write /i, z/ 
and we assume that (4.1) is fulfilled. 

Lemma 4.2. For any 7 > 0, there exists a constant C such that 

i^[-x,x]<C{l + x^). 

Proof. For any positive x, we have 



{l + x)f-Jo (l + y)T 
which, in view of (2.9), is finite. □ 

Observe that the function -0 = ipa,i3 can be written in the form 

1p{x) =l{jfs{x) -iJa{x), 

where 

ippix) =ip^{x) -^pp{x) 



h>0 



; dfi{b,a)- u\ 

\ a a \ Jb<o V a 



a 



dfi{b, a). 



Lemma 4.3. The function tp satisfies (3.8) and (3.9). 
Proof. First we are going to prove 

POO 

(4.3) / \'ipa{x)x\dx < 00. 

Jo 

Choose two small positive numbers di and 82 such that < 61 < 62 < 5, 
where 6 is as in the hypothesis of Theorem 2.2 and take 7 = We may 

assume 7 < 1 . We have 

/ xip+ (x) dx<C e^^ {x) dx 
Jo Jo 
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c 



< c 



-1 



b>OJl 



b>0 
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at /a 

(at-b)/a 
{as+b)/a 

l<as/a J as/a 



dv{s) dt dfi{b, a) 



t^-^dtdu{s) 



{as+b)/a 



t^-Utdu{s) 



I as / a<l<{as+b) / a Jl 

Let us denote the integrals above by / and //, respectively. We have 



a). 



KC 



fe>0 Js>0 



{{as + h)^ - {as)^) du{s) dfi{b, a). 



Observe that for every two positive numbers p, q and every number e with 
< 7 < e < 1 , there exists a positive constant C such that 

{p + qy -p^ <Cp^-'q'. 
Applying this inequality, we may dominate the expression above by 



C 



a 



7-<52/2|„|7-<52/2|L|<52/2 



SxR+ JR+ 

Finally, by the Schwarz inequality, we obtain 



\br''du{s)dn{b,a). 



KC 



-5i/2 



dv{s) 



xR+ 



a-^i (i/i(6,a) 



1/2 



xR+ 



\b\'-d^i{b,a] 



1/2 



and by virtue of (2.10), (2.11), (2.9) and (4.1), the value above is finite. To 
estimate the second integral, we write 

f>(as+6)/o 



II = c 
= c' 
<c' 

<C" 



fe>0 J as / a<l<{as+b) / a Jl 

f I 

6>0 J{a-b)/a<s<a/a 



t'^ ^ dtdv{s) d^t{b,a) 



as 



+ bV 



a 



1 I diy{s) dfi{b, a) 



b>0 J{a-b)/a<s<a/a 

a — b a' 



1 + - 

a 



6\T 



1 diy{s) di2{b, a) 



b^i 



b>0 



a a 



dfi{b, a), 



the last integral being finite by (2.10), (2.11) and Lemma 4.2. Similar ar- 
gument can be used to estimate the integral of ip~, which proves (4.3), and 
moreover, a small modification of the calculations above gives (3.9). 

To prove \xip{x) \ dx < oo, we use the fact that tp = fi*^ f — f . By 
(4.1), we have 



|2;V'(a;)| dx ■ 



^ llogtl 1 

^-^u{at;(3t]dt<C 
t 



at<s<l3t 



dv{s) dt 
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<C / --—dtdv{s)<C s-^du{s)< 

Jo Js/P t'-+'y Jo 



oo. 



Similarly, 

rO 



1 1 



u[aat] (3at\ dtd^A{o) 

dt di'{s) dfiA{a) 



t^+^ 

f3a ps/{aa) ^ 



Js/(f3a) t^^'^ 
Pa 



* f{x)\dx < C 
<C 



<C I I a^s-^dsdv{s) 
Jr+ Jo 

<C J ^a^(^J^ s-^ du{s) + v{l, I3a\^ di^iA{a) 
so the integral is finite by (2.10), (4.1) and Lemma 4.2. □ 

Lemma 4.4. We have 

CIp= I ^{x)dx = log{P/a)-Dl, 
Jr 

where 



D] 



— ;0 

6>o V a . 



dn{b,a) + / u(o; 
Jb<o V 



dfi^b, a). 



Proof. We write 



{^p{x) -'>Paix))dx 



1 



b>oJo t \J(pt-b)/a<s<pt/a 



dv^s) 



-i 



dv{s) ) dtdfi{b, a) 



b>0 



{at—b)/a<s<at/a 
{as+b)/P r{as+b)/a\ i 

- / )-dtdi^{s) 

s>0 \Jas/P J as/a J t 

{as+b)/a ]^ 

- dt di'{s) 

b/a<s<OJ{as+b)/P t 



dfi{b, a) 



= — log(/3/a) [ i^(—-;0 dfi{b,a). 
Jb>o \ a i 

The second part is established in an analogous way. □ 

Repeating the foregoing calculations, one can prove the following: 
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Lemma 4.5. We have 



where 



r2 



Di 



{x + l)V(x) dx = \og{(3/a)iDl + (1 + log(a/3))Z)^), 



+ 

b>OJs>0 Jb<0 J s>~b/a 



log 



as 



as + b 



di>{s) dfi{b, a) 



b>0 J~b/a<s<0 



log(as + b) di'{s) dfiib, a) 



+ 



b<0 JO<s<~b/a 



log(as) di'{s) dfi{b, a). 



Proof of Theorem 2.2. In view of Lemmas 4.2, 4.3, 4.4 and 4.5, we 
apply Theorem 3.3 to prove 

lim = log{(j/a) 



x—^+oo 



r2 ■ 



Analogously, we may describe behavior of the measure v on the negative 
half- line. Namely, the function x i— > /?e^', — ae^] satisfies an appropriate 
Poisson equation and, reasoning as previously, one can prove that all the 
hypotheses of Theorem 3.3 are satisfied. Therefore, we obtain 



x— >+oo 



X 



where 



b>Q 



v[ — ;0 



d^{b, a) 



b<0 



0; 



a). 



Notice -DI = —D\. However, these two constants should be nonnegative, so 
= DI_ = 0. Hence, applying again Theorem 3.3, we prove (2.13) with 



2Di 



and C- 



2D' 



for D\ defined in Lemma 4.5 and 



+ 

6<0"'s<0 Jb>0 Js<-b/a 



log 



as 



as + b 



di>{s) d^{b, a) 



b<0 JO<s<-b/a 



log(|as + b\) dv{s) dfi{b, a) 



+ 



b>0 J-b/a<s<0 



log(|as|) di^{s) dfi{b, a). 



□ 



STOCHASTIC RECURSIONS IN CRITICAL CASE 



17 



Finally, let us observe that the sum of C+ and C_ can be expressed by a 
quite simple formula, namely, assuming (4.1), one can prove 



C+ + C- 







1 log 


as + b 










R 


as 



dvi^s) dn{b, a). 



5. Model due to Letac and extremal random process. Goldie [11] has 
noticed that in the contractive case (i.e., when ElogA < 0) the proof of 
the asymptotic behavior of the invariant measure for the random difference 
equation can be written in a much more general setting and that the same 
ideas can be used for other stochastic recursions (see also Grey [13]). The key 
observation is that the random transformation t ^ At + B can be replaced 
by some other transformations. Although our arguments use the language 
of groups, under more restrictive hypotheses, one can obtain analogous re- 
sults for some other recursions for which the invariant measure is no longer 
probabilistic. 

5.1. Main theorem. In this section we are going to consider the following 
process on M, introduced by Letac [21]: 

Xo = 0, 

(5.1) 

Xn = Bn + An max{C„, Xn-l}, 

where {An, Bn,Cn) are i.i.d. random variables with values in M"*" x M x M"*", 
distributed according to a given measure ^. To simplify the notation, we 
write 

Xn = ^{An,Bn, Cn) o Xn-l- 

This model has been investigated only in the contractive case, Elogyl < 
0. Then there exists a unique stationary probability measure that is, a 
measure satisfying 

ptoi/(/) = l/(/), 

for any positive measurable function /, where 

^^°^{f)= I / f{^{a,h,c)ox)dv{x)d^i{a,h,c). 



Under some further assumptions, the tail of v was described by Goldie [11]. 
We are going to study the critical case 

E[log^] =0. 

As before, we define the measure being the projection of /x onto the first 
coordinate: ha = 7ryi(/i). 

Our main result concerning the Markov chain {^n} is the following. 
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Theorem 5.1. Assume that there exists a positive constant 6 such that: 



(5.2) ElogA = 0; 

(5.3) A^l; 

(5.4) ¥[^{A,B,C) ox = x]<l forallxeR; 

(5.5) HA is spread-out; 

(5.6) KA^<oo and EA~^ < oo; 

(5.7) EjB|^<oo and E|C|^ < oo; 

(5.8) B>6 a.e. 



Then there exists a unique (up to a constant factor) invariant measure v of 
the process {X^}, and moreover, for every positive a < (3, 

lim z/((ax; /?x]) = log(/?/a) • C+, 

a;— >+oo 

where C+ is a positive constant given by the following formula: 
(5.9) C+ = 4/ ^log(^^^^^^)du{s)d^^{a,b,c), 

O" JR+xRxR+JR \ as / 

for 0-2 =E[log2 A]. 

Remark 5.2. The method we use to prove Theorem 5.1 is quite general 
and it is not hard to see that also other examples of stochastic recursions 
can be treated similarly. For example, a special case of the Letac's model, 
which we call extremal random process, defined by the formula 

X'^ = max{AnX'^^i,Dn}, 

where (A„,Z)„) are i.i.d. random variables with values in M"*" x R"*", also 
possesses an invariant Radon measure in the critical case. Under hypothe- 
ses analogous to (5.2)-(5.7) and assuming D > 6 a.s., one can describe the 
behavior at infinity of the invariant measure of {-^^^l- 

5.2. Existence and uniqueness of an invariant measure. To prove exis- 
tence and uniqueness of an invariant measure of the process we 
apply results of Benda [3], who, using ideas of Babillot, Bougerol and Elie 
[1], has investigated locally contractive stochastic dynamical systems. Benda 
has proved that if a stochastic dynamical system {Y^} on M (in fact, he 
worked in much more general settings) satisfies the following conditions: 

• recurrence: for some y, the set of accumulation points of {1^(0;)} is 
nonempty for almost every trajectory lv; 
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• contraction: for every compact set K and every couple of starting points 
liin lK{Yy{u:))\Yy{u:)-Y^{u,)\ = Q, 

n— »+oo 

for almost any trajectory w, 

then there exists a unique (up to a constant) invariant Radon measure u of 
the process {Yn}. 

Proposition 5.3. Under the hypotheses of Theorem 5.1, there exists a 
unique invariant Radon measure v of the process {Xn}. 

Proof. In view of Benda's result, it is enough to justify that the Markov 
chain possesses recurrence and contraction properties. 

Define two autoregressive processes 



and 

= AnXf^_i + Bn + AnCm 



Then both processes satisfy both recurrence and contraction condition [1, 4]. 
Notice 



(5.10) Xl<Xl<X^, 



where X^ is the process defined as in (5.1), but starts from x instead of 0. 
The inequality above immediately implies recurrence of {Xn}- Moreover, for 
any x, ?/ G M, 

\X^n-Xl\<M---An\x-y\ = \Xl-Xl\, 

hence, for compact sets of the form K = [0, M], for every positive constant 
M and almost every trajectory 

\X^M-Xl{u^)\ . 1k{X^M) < \X:iu;)-Xliu;)\ ■ 1,^ (^^(u;)) - 0, 

as n — > oo, which yields the contraction property of {Xn}. □ 
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5.3. Some properties of v. The main result of this section is the fohow- 
ing: 

Proposition 5.4. For any 7 > 0, 

/ : — ;— dvix) < 00. 

7r1 + |x|t 

In fact, this is the only step in the proof where (5.8) is really needed. In 
the proof of the analogous result for the random difference equation (2.9), 
the structure of the group G has been heavily used and the argument cannot 
be applied in our situation. 

For our purpose, we need an explicit formula for z> (cf. [1]). Define a 
random walk on M: 

So = 0, 

(5.11) 

Sn = log{Ai---An), n>l, 
and consider the downward ladder times of Sn- 
Lo = 0, 

(5.12) 

L„ = mf{A; > Ln-i]Sk < Sl„_A- 
Let L = Li. It is known that 
(5.13) -oo<E5l<0 
(see Feller [10]). Next, consider the Markov chain Wn = Xl^. 

Lemma 5.5. There exists a unique invariant probability measure vl of 
the process {Wn}- 

Proof. Observe that the Markov chain {T4^„,} satisfies the following 
stochastic recursion: 

(5.14) 

Wn = max{Z„, MnWn-l + Qn}, 

where (M„, Z„, Qn) are i.i.d. random variables valued in x M x and, 
moreover. Mi =ci e"^^, Zi =^ Xl, Qi =d Xl, where =d denotes equality of 
the corresponding distributions. Therefore, applying (5.13), (5.10) and re- 
sults of Grincevicius [14] and Elie [9] , we obtain 

-oo<E[logMi] <0, 

(5.15) E[log+|Zi|] <oo, 

E[log+|Qi|]<oo. 
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Notice that to obtain the Letac's model, we may write the recursion (5.14) 
in a shghtly different way: 



where 



Wn = Qn + Mn max{W^„_i , Z^} , 

Qn 



zL 



Mn 



Then, under assumption (5.15), existence and uniqueness of an invariant 
probabihty measure was proved by Letac [21] and Goldie [11]. □ 

The following lemma can be deduced from [20]. However, we need it in 
much weaker form than the result proved there and the proof can be con- 
siderable simplified, therefore, we give all the details for the reader's conve- 
nience. 

Lemma 5.6. Let Xi he a sequence of i.i.d. real valued random variables 
such that EXj = 0. Put Sn = ^^=iXi and define the stopping time T = 
min{n : Sn> 0}, then 



E 



T-l 
.n=l 



< oo 



for any 7 > 0. 

Proof. We have 

Define 



E 



T-l 
Ln=l 



:^E[e^'^^r> A;]. 

k=0 



(l)^{s) = Y,s^E[e^^^;T>k], 
k=0 

then, by Spitzer ([26], p. 181, P5 (a), (c)), we have 

1 



(/)^(s) 



' fe=0 



where T' = min{n : S„ > 0}. 

Let fj,' be the distribution of St'- Let Yi,...,Yn be a sequence of i.i.d. 
random variables distributed according to fi' . By [10], 



-00 < EYi < 0. 
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Define S!^ = ^27=1 notice 

fc=0 fc=0 °° 

where U~ =X](a*')*" the Green kernel. The renewal theorem [10] implies 
that, for negative x, U~ {{x,0)) increases linearly, therefore, the expression 
above is finite. □ 

Lemma 5.7. For any positive measurable function f , 



E 



.n=0 



UL{dx) 



Moreover, 



suppz^ C [5, oo). 



Proof. Denote the right-hand side of the foregoing equation by vi. It 
is enough to prove that vi is a Radon measure and it is invariant of the pro- 
cess {Xn}- Invariance can be proved as in [1], page 482. Next, observe that, 
by (5.8), $(j4, C) o [5, oo) c [5, oo) a.s., which implies suppz^^ C [6, oo) and 
suppz^i C [6, oo). To justify that ui is a Radon measure, notice that, for every 
positive constant M, we have 



^^i(l[o,M])=C / E 



X-l 
.n=0 



VL{dx) 



<C J E 



x-i 



J2 l[o,M](e^"^) 



.n=0 



Vhidx) 



< CE 

CM 
< 



L-l 

El[o,M/5](e^") 

.n=0 



•E 



x-i 

.n.=0 



The above value is finite by Lemma 5.6. Finally, because of the uniqueness 
of v, we obtain v = Cvi for some positive constant C . □ 



Proof of Proposition 5.4. Because of (5.10) and (5.8) 

An - ■ ■ AiX + An - ■ ■ A2B1 + ■ ■ ■ + Br, 

>Ai---A 



(5.16) 



\K\>\K\ 
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for every positive x. Therefore, we conclude 



E 



1 



roo 

<C I E 



Ln=0 
L-l 



E 

.n=0 



{i + \x-\p 
1 



duL{x) 
dvL{x) 



and the proposition fohows from Lemma 5.6. □ 



5.4. Proof of Theorem 5.1. We shall proceed as in the proof of Theorem 
2.2. Fix two positive constants a and [3 and define 

/(x) = K(«e^/3e^). 

Then the function / satisfies the following Poisson equation: 

Ji*R f{x) = f{x) + ^{x), 

where 



V'(x) = / / [l(„e-,/3e-](as) 

JR+xRxR+JR 

— l(^ae^^pe'xT^{h + amax{c, s})] dz^(s) (i^(a, 6, c). 
Lemma 5.8. The function if) satisfies assumptions (3.8) and (3.9). 
Proof. Take 7 = (5/4, for (5 as in Theorem 5.1. We have 

e'^'^\ip{x)\dx 

POO P 

i'^^^|l(at,/3t](«s) 

— l(at,/3t] (p + a max(c, s)) \ dv^s) dt dfi{a, b, c) 

i'^"^|l(at,/3t](«s) 

- l(Qt,/3t] (^ + «s) I di^is) dt dfi{a, b, c) 



>cJl 



+ 



-xRxR+ Jc>s Jl 



"^|l(at,/3t](as) 



- l(at,/3t] {b + ac) I di/(s) dt dfi{a, b, c) . 

The fist integral can be estimated as in the proof of Lemma 4.3. The second 
one we dominate by 

o 

*^"Hl(at,/3t] {b + ac) + l(at,f3t] (««)) dt dv{s) dn{a, b, c) 



-xRxK+ Jc>sJl 
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< I / / t^^Ut 

t+xKxIR+ Jc>s \ J{b+ac)/P<t<{b+ac)/a 



t'^ ^ dtj dv[s) dfM{a, b, c) 

as/l3<t<as/a / 



<C I I {{b + acy + {asy)di^{s)dfi{a,b,c) 

-xMxR+ Jc>s 



<c 



<C / s-^'diy{s) ■ 



> + acy ( + («c)^ (^^^ du{s) dfi{a, b, c) 



c''' + a'^c ^) dfi{a, b, c). 



Finally, we use the Schwarz inequality and conclude finiteness of the integral. 
The remaining part of the Lemma can be proved as in lemma 4.3. □ 



Lemma 5.9. We have 



ij){x) dx = 0. 



Proof. Applying the Fubini theorem, we have 

'iIj{x) dx 



b>5 Js>5 



^ ■ l{«i,/3t](as)di 



'i-{at,pt]{'^{a,b,c) o s)dt 



dv{s) diJ,(a, b, c) 



b>5 Js>5 



■dt 



as/(3<t<as/a 



■dt 



(*(a,fe,c)os) /l3<t<{'i>{a,b,c)os) /a 



dv{s) dfi{a, b, c) = 0. 



□ 



Lemma 5.10. We have 

xip{x)dx = log{P/ a) ■ I log 



as 



<I>(a, b,c) o s 



di>{s) dn{a, b, c). 



Proof of Theorem 5.1. In view of previous lemmas, the result follows 
immediately from Theorem 3.3. □ 
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APPENDIX: PROOF OF THEOREM 3.2 

Lemma A.l. The potential A, given by formula (3.2), can he defined 
for every function ij) satisfying assumptions of Theorem 3.2. Moreover, the 
function Aip is continuous. 

Proof. Notice first that 112 * ip is a continuous function, because ip 
is bounded and continuous. Next, a * ip{x) is finite for every x because, 
by (3.17), 

\a*M.'(Pix)\< / \^|J{y)\\a{x -y)\dy 
JM. 

(A.l) <c( [ \ip{y)\dy + \x\ f my)\dy+ f \tl;{y)y\dy 

\Jm. JR JR 

<C(V')(l + |x|) <CX). 

Finally, to prove continuity of the function a*ij), fix x G R and consider a 
sequence {xn} tending to x. Put h,n{y) = ^p{y)a{xn — y)- Then by (3.17), all 
the functions are dominated by h{y) = \tp{y)\ ■ {\a{x — y)\ + C), for an 
appropriate large constant C, which is an integrable function. Therefore, by 
the Lebesgue theorem and using continuity of the function a, 

lim a*M^(x„)= lim / hn{y)dy 
= / ip{y)a{x - y) dy 

JR 

= a*R^/^(x). □ 

Lemma A. 11. For every function tp satisfying assumptions of Theorem 
3.2 and x S M, the following Poisson equation is fulfilled: 

Jl *M Ail:{x) = A'4){x) + ii{x). 

Proof. The foregoing Poisson equation is satisfied when ip is an ele- 
ment of $ (Port and Stone [24], Theorem 10.1). Without any loss of gener- 
ality, we may assume ip'>^. Then take any sequence ipn of positive contin- 
uous, compactly supported functions, tending pointwise to ij) and satisfying 
V'n(x) < il^{x) for every x and n. Then 

p. *K Atpn{x) = Atpnix) + tpn{x). 

Notice that Aipn tends pointwise to Aip, which is just a consequence of the 
Lebesgue theorem. Finally, observe, that by (A.l), Aip(x) can be bounded 
by C{\x\ + 1) for some constant C, therefore, /2 * Aip is well defined and 
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applying again the Lebesgue theorem, we conclude that /2 * Aipn tends to 
p,*A'ip. □ 

Let us fix a function /, which is a solution of the Poisson equation (3.1), 
for some function ip satisfying hypotheses of Theorem 3.2. Suppose g is 
a continuous, compactly supported function, such that J{g) = 1. Define a 
function 

(A.2) h{x)=f{x) + J{'ip)Ag{x)-A'ip{x), x eR. 

Lemma A. 12. If the function f is bounded from below, then also the 
function h is bounded from below. 

Proof. It is enough to show that J{'ip)Ag — Aip is a bounded function. 
By Lemma A.l, it is a continuous function, therefore, it is enough to justify 
that the limits 

lim (a*R {J{ilj)g -'4)){x)) 

X— »±oo 

exist and are finite. We have, by (3.17), 
lim (a *K {J{'>p)g - ip){x)) 



= lim {a{x-y)-a{x)){J{'il))g{y)-iij{y))dy. 

Observe that, because of (3.17), 

\a{x-y)-a{x)\\J{il})g{y) - i}{y)\ 

<C(|y| + i)|J(V)<7(y)-V'(y)l, 

which is an integrable function. Therefore, by the Lebesgue theorem, we 
obtain 

lim (a*R(J(V')5f-V')(2;)) 

X— » + oo 

= -o-"^/ y{JWgiy) -i^iy))dy <oo, 

Jr 

which proves the lemma. □ 

Proof of Theorem 3.2. Define h as in (A.2), then h is bounded from 
below (Lemma A. 12) and by Lemma A. 11 satisfies 

li*^h{x) = h{x) + J{il))g{x), X e M. 

Therefore, by Theorem 3.1, 

cJhp) 

h{x) = J Ag{x) -\ 2 — x + d, dx a.s. 
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Hence 



f{x) = h{x) - J{i))Ag{x) + Ail^ix) 



(A.3) 



= A^{x) + ^-^ + d, 



dx a.s. 



□ 
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